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■ In this work we show that, under given conditions, parametric Backhand trans- 

CN . formations (BTs) for a hnite dimensional integrable system can be interpreted as 

the solution to the equations of motion defined by an associated non- autonomous 
Hamiltonian. The two systems share the same constants of motion. This obser- 
£Sj . vation lead to the identification of the Hamiltonian interpolating the iteration 

of the discrete map defined by the transformations, that indeed will be a linear 
combination of the integrals appearing in the spectral curve of the Lax matrix. 
^\ ■ An application to the Toda periodic lattice is given. 
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p 1 Introduction 

^ ■ In the last years there has been a huge interest in the theory of integrable discretization. 

With "integrable discretization" usually is meant (see e.g. [16]) a discrete analogue of 
an integrable system of equations (H being the Hamiltonian, possibly a function of all 
', the N integrals of the model and {•, •} being the Poisson brackets): 

z i = fi(z)={H,z i } (1) 



that is a map (implicit or explicit) 



depending on some parameter \x such that approximates, at least at first order, the 
flow [[] (i.e. Zi = Zi + ^fi(z) + 0(/i 2 )); that preserves, in the limit \i — > 0, the Poisson 
structure {■, ■} and such that possesses N conserved quantities independent and in in- 
volution, possibly depending on /i but coinciding with the continuous ones in the limit 
li ->■ 0. 

A great number of results, based on various methods, appeared in the literature in the 
last two decades on discretization of finite dimensional systems fitting the previous def- 
inition pQ [3] [2] [5] [TQ] [H] [12] [13] (see also [16] and references therein). Among the 
various techniques, the Backlund transformations approach appear to be noteworthy: 
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indeed, despite the fact to be "non-algorithmic" (unlike other methods, e.g. the "split- 
ting method"), it possesses characteristics, such as explicitness, exact preservation of 
all the integrals and canonicity (see e.g. [7]) deserving special consideration. 
The aim of this paper is to show how Backlund transformations for finite dimensional 
systems represent, provided the "spectrality property" [5] holds, an exact discretization 
of the underlying continuous integrable model. With the adjective exact we mean the 
equivalence of the trajectories, at all times and at all orders, of the discrete flow defined 
by the transformations and the trajectories of the continuous flow of the model. 
In Section (j2j) we consider an integrable system described by a 2x2 Lax matrix pos- 
sessing a set of (parametric) Backlund transformation with the spectrality property. 
We will show how these transformations can be seen as the integral curves of a system 
of non-autonomous Hamiltonian equations sharing the same conserved quantities with 
the ancestor model. Then it will be pointed out that indeed the recurrences defined 
by the transformations lie on the trajectories of the continuous system defined by the 
Hamiltonian 

n = J2 c kH k (2) 

k 

where are the conserved quantities of the continuous model. The constants Cf. will 
depend on the free parameters appearing in the transformations. At the end of the 
Section we will discuss the "usefulness" of the formulae obtained. 
In Section Q we will give an application of what stated in Section (J2J) with a two- 
parameters BTs on the periodic Toda lattice. The last paragraph will present some 
numerical examples on this model. 

2 BTs & integrable discretization: an Hamiltonian 
approach 

In this Section we consider an integrable system with a 2x2 Lax matrix L(X), A being 
the spectral parameter. The conserved quantities (independent and in involution) of 
the model can be obtained by the coefficients of some expansion in A of the determinant 
and of the trace of L(\). So the characteristic equation: 

det(L(A) - vt) = v 2 - Tr(L(X))v + det(L(A)) = 

defines the eigenvalues of L(X), v\ and i>2, in terms of the conserved quantities and A. 
It is convenient, but not necessary, to think about our model as written in terms of 
canonical variables, say {pi,qi}f =1 , so that we have the usual Poisson brackets: 

{pi, qj} = 5ij {pi,Pj} = {qi, qj} = 

and the involutivity conditions on the conserved quantities Hf. 

{H i: Hj} = Vi,j = l...N (3) 



2 



A set of Backhand transformations for the system is simply a set of canonical transfor- 
mations 

Pi=>Pi = fi(p,q) () 

q l ^q i = 9i{p,q) 

from (pi, qi)f =l to a new set (pi, qi)f =1 such that iV functions (Hi)f =l are invariant under 
the transformations, that is Hi(p,q) = Hi(p,q) (here and after p and q stand for the 
sets (pi)f =l and (qi)f =1 respectively). From this definition indeed the classical property 
of BTs (or better auto-BTs) to send solutions of the equations of motion into solutions 
is easily proved, since the two sets of canonically conjugate variables satisfy the same 
system of equations. The application of BTs in numerical analysis comes from the 
observation that actually the relations (jlj) define explicit recurrences just by thinking 
at tilded variables as the untilded ones but computed at the next time step (in the 
following notation pi is equivalent to (pi)o, Pi is equivalent to (pi)\ and so on): 

: fi((p)n, (q)n) ^ 
9i((p)n, (q)n) 

As pointed out by Kuznetsov and Sklyanin [8], the invariance of the conserved quanti- 
ties, and then of the spectrum of L(X), implies the existence of a dressing or Darboux 
matrix D(X) intertwining the two matrices, L(X) and L(X) = L(X,p,q): 

L{X)D(X) = D(X)L{X) 

Also, because the transformations are canonical, it exists a generating function, say 
F(q, q), such that 

dF(q,q) 




Pi 



Pi 



dqi 

dF(q, q) 



Of particular interest are the BTs depending on a parameter, say li. In this case one 
can always think that the transformations are connected to the identity, in the sense 
that for a value of the parameter, that can be always be chosen to be 0, one has p = p 
and q = q. Indeed if the original transformations are not connected to the identity we 
can obtain a new set of transformations that, by construction, are the identity trans- 
formations for li = 0, by the following very simple argument (see also [7]). 
Let the dressing matrix defining the transformations be D(X, X\) where Ai is the param- 
eter of the transformations. Since BTs are canonical transformations we can compose 
a transformation with parameter Ai and a transformation with parameter A2 to obtain 
a two-parameters set of BTs. So composing the transformation with Ai = Ao + [i with 
its inverse transformation but calculated at A2 = Ao — \i will give a transformation that 
in the limit li — > goes to the identity. The only question is about the concrete invert- 
ibility of the map, but having the expression of D(X, X\) the inverse transformation is 
found just by the adjucate of the matrix D(X, Ai), say D(X, X\). The dressing matrix 
for the composed transformation then will be D(X, Aq,aO = D(X, A2) • D(X,X\). We 
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will give an example of such construction in Section (J3J) 

It is quite natural to interpret the parameter fx as a sort of "time" with respect to which 
the BTs evolve. As a matter of fact Backlund transformations are parametric canonical 
transformations; this means that they are the integral curves (whose parameter is fx) of 
an Hamiltonian vector field since, at least locally, there exist an Hamiltonian K(p, q, fx) 
such that (see for example [I] for the proof): 

dK 




dpi (6) 
dqi 

As now we will show, if the so-called spectrality property, introduced by Kuznetsov 
and Sklyanin [8] holds, then it is possible to find explicitly the function K(p,q,fx); 
this function is indeed related to the eigenvalues of the Lax matrix and so contains all 
the conserved quantities of the system, but it depends also on the parameter(s) of the 
transformations: BTs are then the integral curves of a non-autonomous Hamiltonian 
vector field (the time being fx). 

In order to clarify the previous discussion let us recall what the spectrality property 
is. Suppose to have a set of BTs depending on a parameter fx: 

Pi = fi(p,q,v) 
qi = 9i(p,q,fi>) 

Assuming the inverse function theorem applies, we can say that there exists a generating 
function F such that: 

dF(q,q,fx) 



Pi 



Pi 



dqi 
dF(q,q,fx) 



The spectrality property says that (a function of) the canonically conjugate variable 
with respect to fx (say $) and fx itself lie on the spectral curve: 



dfx 



=^ det(L(//)-/($)U) = (7) 

<?=<I(p,<?,At) 



It seems that this property is "universal" in the sense of being shared by a large class 
of models (Toda [15], Ruijsenaars [12], Henon-Heiles [6], Gaudin [7], Kirchhoff top [H], 
Mumford systems [9]...). 

According to some authors it could be better understood by taking an algebraic, 
or geometric, point of view. Here we want just to give a mechanical (Hamiltonian) 
interpretation of this property. Our observation is very simple: the "variable" $ is 
canonically conjugated to fx; furthermore, when expressed in terms of the variables p 



^ee e.g. |15j where the author writes "the property .. .needs more research to uncover its algebraic 
and geometric meaning" . 
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and q, it is a function of the conserved quantities because lies on the spectral curve. 
Enlarging the phase space by adding the couple ($, /i) we can write dF as the difference 
between two Poincare Cartan forms: 



dF =} pidqi - pidqt + 
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This means that if (p, q) are the integral curves generated by the Hamiltonian T(p, q, /i), 
then p and q are the integral curves generated by the Hamiltonian 



where T is the transform of T(p, q, /i) when expressed in terms of p and q. Notice that 
the function $, thanks to spectrality, is invariant under the transformation, so $ = $. 
Now let us make this further observation: as explained some lines above, it is always 
possible to have a BTs connected to the identity, so that: 



Since we have the freedom to choose the function T(p,q,fl), we can choose it to be 
zero. With this choice (pi,qi)f = i are just constants (the initial conditions). From (Q 
and f lTUj) then readily follows the following 

Theorem 1 Assuming the spectrality property ([?]) holds, the BTs (pi,qi)iLi ar ^ the 
integral curves of the following system of (non- autonomous) Hamiltonian equations: 



corresponding to the initial conditions [TO]) . 

By this point of view obtaining a set of explicit BTs transformations with the spectral- 
ity property for an integrable model means to have the explicit general solution of a 
non-autonomous Hamiltonian system sharing the same conserved quantities with the 
model. So we are considering two set of flows: the first set is given by all the flows 
defined by the integrable model, the second set is given by the flow defined by $. In 
the following we aim to clarify the relations among these flows in order to point out a 
concrete application of Theorem [TJ 

The key point is that the the flow defined by (TTTT) possesses the same conserved quanti- 
ties ([3]) as the corresponding integrable model defined by the Lax matrix L(X). Being 
the model integrable (by definition, in the sense of Liouville), there exist action-angle 
variables (Ik,Qk)k=n such that the action variables Ik are functions only of the con- 
served quantities and the couples (Ik, Ok) are canonically conjugate [3]: 



K{p,q,y) = T(p(p,q,fi),q(p, q,fJ,),li) + $(p,q,(j) 



(9) 



(10) 




(11) 



{Ik, 0j} = 5kj 



{Ik,Ij} = 0, 



{O k ,0j} = O 



(12) 



5 



Actually we do not need that all the variables 9j are really "angles", it could be also 
a splitting among bounded and unbounded variables (say 9k G T m ,k = l..,m and 
9k G R^"" 1 , k = m + 1...N). The key request is the canonicity (fT2|) (obviously together 
with the conservation of each of Ij). Having this in mind, for the sake of simplicity we 
continue to use indistinctly the term "angle" in what follows. 

Let us see what happens to the flow defined by $ fill I) in these coordinates. When 
expressed in terms of the set (I k , 9k)k =1 the Hamiltonian $(p, q, fi) = $(j>, q, /i) reduces 
to a function of action coordinates and only: 

$(p,g,/x) = $(J,aO = Hi,!*) 
The system (fTTl) then becomes: 

(di 3 _ a$(/, / u)_ n 
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dfi d9 

d§ j _d${f,IA) _ Q , d 

I dfi dl 



j U1 i Jo 



d - ~ 

oh Jo 



(13) 



In these coordinates it is clear what happens to the recurrences (J5]) defined by the BTs. 
Their explicit solution is indeed given by: 

(Ij) n = (ij)o 

d - ~ (14) 
(9 j ) n = (9 1 ) + n- r / 1 ) 

olj Jo 

Replacing n by t we see that equations ([Till represent the exact time discretization of 
the continuous system governed by the Hamiltonian: 

H= r$d[i (15) 
Jo 

Notice that there is no need to find explicitly the action-angle coordinates: the result 
(1151) is independent of which coordinates we are expressing the equations of motion; 
we could also forget at its derivation and state simply the following 

Theorem 2 If the spectrality property holds, the discrete trajectories defined by the 
iteration of the recurrences 

(Pi) n +i = fi(pn,q n ,n) 
{Pi)n+\ = gi(p n ,q n ,n) 

given by the BTs coincide with the trajectories of the Hamiltonian system defined by 
the Hamiltonian: 



where <3> is defined by the spectrality property |?p. 
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So, by one hand the BTs are the integral curves of the system of Hamiltonian equa- 
tions generated by $, by the other hand they are the exact time discretization (i.e. 
preserving the trajectories) of the Hamiltonian equations generated by % = J <&dfx. 
When expressed in canonical coordinates {pi, qi)f =1 , $ is a function of the conserved 
quantities Hi of the model: $ = Q(H,fx). The BTs are then the discretization of 
the flows corresponding to linear combinations of these conserved quantities; indeed 
for any function $(p, q) on the phase space we can write: 

d = {% 3} = c dH k , 3}, where c k = -|- T ®dfx (16) 

k ^ 

Let us make some remarks about the "usefulness" of the previous formulae. Having 
in mind a physical model, usually one looks only at a discretization in some particular 
direction, corresponding to a given choice of the constants Ck appearing in ( fill . So the 
problem is about the freedom we have to choose the values of these constants. 
As noted before, we always have the possibility to compose BTs writing products of 
dressing matrices with different parameters. The result will be a BTs depending on a 
set of free parameters. Also, the generating function F © and the spectral function $ 
(and then the constants Ck themselves) will depend on these parameters. If there are 
enough parameters (that is N, the number of degrees of freedom of the system), the 
problem is then to invert the integrals for the constants Ck- 

< *=^r** <i7) 

For example for algebraic complete integrable systems one expects these integrals be 
hyperelliptic integrals. Notice however that for any fixed orbit and for a given set of 
the constants c k , the integrals could be inverted numerically. An example will be given 
in the paragraph (13. 2p . 

The parameter fx plays a role of a time not only for what concerns Theorem (CQ), but 
as it is clear also from (I14p . at smaller values of fx there corresponds smaller values in 
the time step discretization. 

In the next Section we are going to present an application of the previous formulae. 



3 An example from the Toda lattice 

Let us consider the BTs of the periodic Toda lattice. The Lax matrix of the model is 
given by L(A) = In(X) ■ ■ ■ h(X) where the local matrices lj(X) are (see e.g. 



«*> - ( A e-,f T ) us) 

The variables (pj,qj)f=i are couples of canonically conjugate variables. Also, for later 
use, it is useful to define the four functions (A(X),B(X),C(X),D(X)) as the elements 
of the lax matrix L(X): 

~ ( $ ) M> 
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The Poisson brackets among any combinations of these functions are defined by the 
following r-matrix structure: 



{L(X),L((i)} = [r(X - (i),L{\) ® L(/j,)] , (20) 
where r(A) is explicitly written as: 



r (A) = -i 



/ 1 \ 

10 

10 

\ 1 ) 



(21) 





-A(/i)B(A) 


A 


- H 


C(X)B{fi) 


- C{n)B{X) 


A 


-n 


D(A)S(/i) 


- D(n)B(X) 



The corresponding Poisson brackets among the functions (A(X), B(X), C(A), D(X)) read 
as: 

{A(X),A(»)} = (fl(A),BW) = {C(A),CM} = {C(A),DM} = (22) 

{A W ,B M } = — ; M(A),co,)} = ^ c ' A »-^ A » c "" 

A — /! 

(23) 

M(A) , so,,} = :~ ; {b( a), = ^j:^ 1 

(24) 

A — fl A — /i 

(25) 

Note that the determinant of each local Lax matrices lj(X) is 1 so that the determinant 
of L(A) = A(X)D(X) — B(X)C(X) is also 1. The hyperelliptic curve defined by the 
characteristic polynomial det(i> — L(X)) = 0, that is 

v 2 - tx{L{X))v + 1 = 

is constant with respect to the Toda flows. All the commuting Hamiltonians Hj of the 
system can be obtained by the coefficients of the expansion of tr(L(A)): 

N 

tr(L(A)) = A(X) + D(X) = H k \ k (26) 

fc=0 

Note that there are N Hamiltonians because of = 1. 

Following [15], we can construct a class of BTs for this system considering a similarity 
transformation through a matrix D(X) (the Darboux matrix) defining a map from the 
variables {{pj,qj)}f=i to the variables {{pjAj)}f=i- 

L(X,p,q)D(X) = D(X)L(X,p,q) (27) 

Due to the splitting L(X) = Zjv(A) ■ • • h(X) of the Lax matrix into local matrices, it is 
possible to recursively define a set of local matrices Dj(X) by the following relation: 

lj(X,p,q)dj{X) = d j+1 (X)lj(X,p, q) (28) 
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so that the transformation (127)1 is fulfilled by the position D(X) = di(X). The simplest of 
the matrices dj corresponds to the choice of linear dependence on the spectral parameter 
H3|: 

d i w = ( x ~ c p j ajPi 7 J ) (29) 

where £ is believed to be a parameter and ccj and j3j depend on the dynamical variables. 
The dependences of ctj and f3j are immediately fixed by the off-diagonal part of (EE}, 
giving: 

' ft+1 T' (30) 

With these positions, the similarity transformation (128]) defines the following implicit 
map: 

Pi (<'' " ■ - " ■ C) 

Although implicit, it is very useful that the map has been written in this form in order 
to prove its canonicity. Indeed it can be easily checked that eqs. ( 13"Tj) can be found by 
the generating function F((, q, q) defined as follows: 

BP flF 1 

Pj=QT> p i = -W- i? = E(Cfe-^) + e"-e^^) (32) 

3 3 f. 

It is not an accident however that the matrices dj (and hence D — d\) have a determi- 
nant vanishing at A = Q. In this case indeed the transformations (13"T|) can be plainly 
inverted (in the sense to write explicit transformations). Indeed D(\) is given by: 

D(X) = ( X ~ C ~ aiPl "f) (33) 

We already know that a.\ = e qi and j3\ = e~^ N . But, since det (D(Q) = 0, D(Q has a 
(unique, up to an overall factor) kernel, given by \Q >= (1, — /?i) T . The combination 
of the unicity of this kernel with (J27l) gives that \Q > is also an eigenvector of L(Q. 
Calling v(Q the corresponding eigenvalue, we have: 

f A(C) - f3 1 B(C) = v(C) 

\ C(C) - A£>(C) = -v(C)Pi 1 J 

After the elimination of v(Q we are left with (3fB(() - Pi(A(C) ~ D(Q) - C(() = 
defining f3\ in terms of only the untilded dynamical variables. At this point the relation 
(127]) gives explicit transformations. 

So we have BTs with a free parameter and, as explained in Section (j2j), we can ask 
about the flow associated with this parameter. We remember that, since the map 
(jHTj) defines canonical transformations, there exist (at least locally) an Hamiltonian 
K(p> <?> 0) sucn that 

dq dK 

% _ %_K_ (35) 
d( dq 
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The flow defined by (}3~Tj) is however not connected to the identity. In order to obtain a 
flow connected to the identity, as explained in Section (j2J), we can compose two BTs, 
the first with parameter £ = Ai and the second given by the inverse transformations of 
the first one but evaluated at £ = A2. It is clear then that when A2 — > Ai the identity 
map is recovered. So let us consider the inverse transformations of (131 p . They are 
defined by the inverse Darboux matrix D~ l but since ( 12 7p is homogeneous, we can take 
the adjugate matrix of D, say D. So we have 

L(X,p,q)D(X) = D(X)L(X,p,q) (36) 

with 

^iM-c'-J (37) 

We know that ol\ = e gi and f3\ = e~ qN but in order to have explicit transformations 
we need to write a\ in terms of tilded variables. Note that -D(C) has a kernel, given by 
\Ct >= (ai, — 1) T that is also an eigenvector of L(£) = L(£,p, q). So we have: 

r i(o«i - b(o = w(o ai 

\c(C)ai-D(Q = -w(Q 

After the elimination of v(() we are left with (ai) 2 C(C) + «i(i(C) - D(()) ~ B(C) = °- 
Note also that since det L = det L = 1 there are two possibilities for the eigenvalue 
iw(C): or w(C) = u(C) or W (C) = t '~ 1 (C)) where ?;(£) is defined by (134)) . In order to have 
the inverse transformations of the previous one the correct choice is w(Q = w~ 1 (C) a s 
now we will show. This will fix the relative sign of the square roots appearing in the 
solutions of the quadratic equations defining /3i and o^. 

The element (1, 1) of ([35]) gives A(() - PiB{£) = A{() + aiC((). But from ([21 we 
have A(0 — fi\B(Q = v(£). Now is simple to see that the equivalence w(Q = v _1 (£), 
or more explicitly: 

(D(Q-a 1 C(t))(A(Q + a 1 C(C)) = l 
is just the equation defining a.\ because of AD = BC + 1. 

So we are ready to compose the two maps. We can define a new Darboux matrix 
T>(\, Ai, A 2 ) given by the product D(X, \ 2 )D(\, Ai): 

LD(X, AO = D(X, Xi)L, ~LD{X, A 2 ) = D(X, A 2 )L, — ► 
— ► LD(X, X 2 )D(X, A x ) = D(A, A 2 )D(A, A X )L 

with: 

^ ^ - ( 4 a - C- xm ) - ( A " A ;r " yi t ) < 4o > 



From the relations ( 1391) one has y\ — y 2 — e~ qN . So we set y\ = y 2 = y. Also the 
product D(X, X 2 )D(X, Ai) depends only on the difference x\ — x 2 , so we set x\ — x 2 = x. 
The result is: 
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From now on we consider directly the transformations given by this Darboux matrix, 
so omitting the intermediate step we look at: 

L(X)V(\, Ax, A 2 ) = V(X, Ax, A 2 )L(A) (42) 

The matrix ( HIT) is obviously degenerate for A = Ai and A = A 2 . The respective kernels 
are given by \Q± >= (1, — y) T and |f2 2 >= A 2 — Ai— xy) T . These are also eigenvectors 
of L(Xi) and L(X 2 ) with eigenvalues v(Xi) and v~ 1 (A 2 ): 

L(X 1 )\n 1 >=v(X 1 )\n 1 >, L(X 2 )\n 2 >=ir\X 2 )\tt 2 > (43) 

We can explicitly write: 

4(Ai) + P(Ai)- 7 (Ai) A(A 2 ) + D(X 2 ) + 7 (A 2 ) 

^(Ai) = , v(X 2 ) - 



2 v ZJ 2 (44) 

with 7(A) 2 = (A(X) + -D(A)) 2 - 4 



From relations (|43|) the variables x and y are then easily found: 
A(X 1 )-D(X 1 )+ 1 (X 1 ) 



X 



2B(X 1 ) 

2B(X 1 )B(X 2 )(X 2 -X 1 



(45) 



B(X 2 )(A(X 1 ) - D(X l ) + 7 (Ai)) - 5(A 1 )(A(A 2 ) - D(A 2 ) - 7 (A 2 )) 
The relations above enable us to obtain in explicit form the BTs. Posing 

Ai = A + /i, A 2 = A - /i (46) 
we see that in the limit fi — > the identity transformation is recovered. 

3.1 Generating functions and Hamiltonian flows 

In order to find the generating function for the composed transformation produced by 
the Darboux matrix (]4ip it is possible to proceed as for the case of one parameter 
transformations. So, let us find the implicit transformations p = p(q, q, /i) and p = 
p(q, q, (i) in order to obtain the generating function F as the solution of the system: 

dF ~ dF 

p > = %' * = ~Wi (47) 

Again we look at the relation among the local Lax matrices: 

^■(A,p,g)^(A) = Vj+iWljiX^q) (48) 
where now the matrices Dj are given by: 

V %(Ai - A 2 + XjUj) X-X 2 + x j y j J 
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Recall that d\ = T> so that x,\ = x and y\ = y. From the off-diagonal elements of ([48]) 
we obtain the expressions of Xj and y,j in terms of the dynamical variables qu and qk 
as: 

y i+ i(2/i + x j+1 y j+1 ) = e~* - e~* (50) 



It is useful to introduce the variable rjj defined by: 

Vj = XjVj + V 

Indeed from: 

x j+1 y j+1 (2n + x j+1 y j+l ) = (e* +1 - e* +1 ) (e~* - e - *) = 

We see that 77? = /x 2 + Wj. With the help of these relations, the remaining equations 
gHD give: 



Pj = -A 
Pi = -A 



After some calculations it is not difficult to check that the generating function F for 
these transformations is given by: 



(^ ln 



Vk + At 



2% - A (g fe - <? fc ) (51) 



Now, as explained in Section |2l we can consider an extension of the phase space defining 
the variable conjugate to /i: 

(52) 

By an Hamiltonian point of view $, when evaluated at q — q(p,q,p), is no other 
than the new Hamiltonian flow defined by the canonical transformations engendered 
by the generating function F. So if the variables (pi, qi)f = i represent the integral curves 
generated by an Hamiltonian T(p, q, //), the variables (pi, qi)i = \ are the integral curves 
generated by the new Hamiltonian: 

K(P, Q, A 4 ) = T(p, q,fj) + $ (p, q, /i) (53) 

Also, taking T = 0, we can consider the values of (Pi-,qi)f = i as the set of our initial 
values; the BTs are then the general solution of the system of equations given by the 
Hamiltonian $. 

The function <3>(j5, q, fi) assumes the same value when evaluated at (p, q) or (p, q) due 
to the equivalence A(X) + D(X) = A(X) + D(X) (recall that A + D is the generating 
function of all the conserved quantities of the system, see eq. So 

$(p, 9) A 4 ) = ?(P, 9, A 1 ); A 4 ) (54) 
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At this point it is convenient to calculate explicitly the function $ in (153]) . 
From (1521) we have that: 



u \ 



k 

In order to write $ explicitly in terms of {pi,qi)f = i we note that the matrices Dj (149|) 
are degenerate when A = A x and A = A 2 . The kernels of Oj(Ai) are given by |o>,,-(Ai) >= 
(xj, \i — r]j) T and the relations give: 

^■(AOIu^Ai) >= ^ + i|w i+1 (Ai) > 

for some function More explicitly: 

(pj + Ai)xj - (/i - 77^) = v} +1 x j+ i (55) 

But from the combination of the elements (2, 1) and (2, 2) of we have: 

( Pj + X) Xj = e q > (A - A - + (A - A + %+i) (56) 

Evaluating the above equation at A = Ai and confronting with (I55p leads to the iden- 
tification of the function isj +l as: 

i _ e*(/i + ?/ i+ i) 



Repeating the same argument for the kernels \uj(X 2 ) >= (—Xj,/i + i]j) T of dj(X 2 ), we 
can write for some function Vj +1 : 

^•(A 2 )|wj(A 2 ) >= Uj +1 \u j+1 (X 2 ) > 

and again with the help of ( |56l) evaluated at A = A 2 we find: 

.2 _ e 9_J (?7i+i - AO 



^+1 



X J+ i 



Now due to the factorization of £(A) in local matrices it is straightforward to show 
that the products Ylk u l an d Ylk v \ are eigenvalues of L{X\) and L(A 2 ) respectively. 
But from the analysis of the eigenvalue of L{X\) and L(X 2 ) of the previous section we 
know that, when expressed in terms of the untilded variables we have (see f H3|) and 

I): 

L(Ai)|fii>=v(Ai)|fii>, L(X 2 )\n 2 >=v- 1 (X 2 )\n 2 > (57) 



Collecting all together we see that ^ as in ( 1521) is given by: 



dF 



In 



= In 

q=q(p,q,^i) \ 1 Ifc "A; 

A(A X ) + D(A X ) - 7 (Ai) A(X 2 ) + D(A 2 ) - 7 (A 2 



TLX 
11^? 



= ln (v(AiMA 2 )) = 
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This is the spectrality property of Kuznetsov and Sklyanin (see e.g. [8]): a function of 
the canonically conjugate variable with respect to the parameter of the transformations 
and the parameter itself lie on the spectral curve: det(f(v) — L(fi)) = 0. The result 
can be also rewritten in terms of inverse hyperbolic function arccosh(-): 



dF 
dpi 



-arccosht^M+giM) _ arccosh( ^M±^M) (59) 



<?=<j(p,<?>a0 



The BTs given by ( l4Tj) are then the integral curves of the flow generated by the Hamil- 
tonian (1591) . 

As an example let us take the simplest case, that is N = 2. 
The trace of the Lax matrix is now given by: 



A(X) + B(X) = A 2 + (p 1 + p 2 )X + P1P2 - 2cosh(g! - q 2 



(60) 



In this case it is also quite straightforward to find explicitly a set of action-angle 
variables. A choice is given by the following relations: 



Pi +P2 



61 = qi + <?2, 
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1-2 = 2cosh(gi — q 2 
F( P2 ~ P1 k) 



Pi -P2 



(61) 



where F(z, k) is the complete elliptic integral of the first kind and the modulus k is 
given by k 2 = Now expression fl60|) is written as A(X) + B(X) = (A + /1) 2 — I 2 
According to Theorem ([2]) the recurrences defined by the BTs §5§ discretize the con- 
tinuous flow corresponding to the Hamiltonian: 



N=2 



so that 







0i = 61 + — — / arccosh 



arccosh 



(A + h) 



dX 



(A + h) 



6, 



02 + ^ 

dh 



Ao-M 



arccosh 



2 J\ +fi 



(X + h) 



dX 
dX 



giving 



6*i = 1 + arccosh 
02 — 02 + 



(A -/i + /i) 2 -/ 2 



arccosh 



(A + n + hf - h 



r( A -M+_J l fe) _ F( Ao + M+_J 1 )A;) 



(62) 



(63) 
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In these coordinates the BTs are linearized. Indeed if we denote, as in (fT4l) . with (Ij) 
and (6j) n the n th iterates of the variables Ij and 6j, we obtain from (163|) : 



f (h)n = (h) 

/ l\ .. I J \2 T \ / I \ , .. , J \2 7 

2 



(A - ii + hf - h\ , /(A + /i + /i) 2 -/ 2 



(#i)n = (#1)0 + narccosh f ) — narccosh 

(64) 

Also, confronting with the definition of 6 2 (IBTj) . it becomes clear that the BTs corre- 
spond to addition formulae for elliptic (or hyper-elliptic in the case N > 2) integrals. 
Indeed we can write: 



F(P, k) = F(P, k) + F(\-,k) - F(X+, k) (65) 

where P = J^j£= , P = f 2 ,7 Pl n and A± = A °^ = r M+ r f 1 . For analogous results on addition 
formulae for WeierstraB p function or Jacobi elliptic functions see respectively [9] and 



3.2 Numerics 

The results of Section (]3]) together with Theorem ((2]) lead to the following result: the 
BTs for the Toda discretize the eqs. of motion given by the following Hamiltonian: 

n = - / A0+ "arccosh( Tr(L(A)) )dA (66) 

where we recall that the trace of L(X) is the generating function of all the integrals: 

N 

tr(L(A)) = A(X) + D(X) = ^ H k X\ H N = 1 (67) 

fc=0 

so that the corresponding continuous eqs. of motion, for any function $ of the phase 
space, is: 

$ = {H, $} = V c k {H k , /}, where c k = - / =dX (68) 

4" An-« A/Tr(L(A))2-4 1 ^ 



Let us take the explicit case N = 3. The three integrals of motion are given by: 

H 2 = V1V2 + V2V2 + PsPi ~ e qi ~ qa - e q2 ~ qi - e 93 " 92 (69) 
H 3 = Pl p 2 p 3 - p 2 e qi ~ qs - p 3 e q2 ~ qi - pie 93-92 

Given the formulae (1681) it is possible to plot the continuous flow corresponding to some 
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Figure 2: The values of [i are 5 (top) and 2 (bottom). The corresponding values for the 
constants (c ,c 1 ,c 2 ) are (-0.005,-0.048,-0.512 and (-0.0018,-0.0183,-0.1849) 
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choice of the constants Ck and the corresponding discrete flow given by the recurrences 
defined by the explicit BTs (from (I42p . (I45p and (14*6]) ). In order to make a comparison 
one has to choose the same initial conditions for the two flows. Once they are fixed 
(so that the conserved quantities H\, H 2 and H 3 are given), it is possible to calculate, 
numerically or explicitly, the values of the constants knowing those of Ao and ji. 
Conversely, if there were three free parameters instead of two in the BTs, it would be 
possible to calculate their values knowing those of the constants c^. Let in our case 
choose the values Ao = 10 and \x = 5. A plot of the continuous flow given by ([151) is 
given in figure (OQ) on the left. On the right there are the first 3500 iterations of the 
BTs with the same initial conditions (both initial conditions and the values of the c^'s 
are reported in the caption). The variables on the axes are (p3, q 2 — Q, qz, — Q), where 
the center of mass Q = qi+q ^ +qz corresponding to a linear motion, has been subtracted 
to reveal the structure of a 2-torus, in accordance with the Liouville- Arnold theorem. 
In figure ([2]) we report an example of the continuous trajectories (in grey) and the dis- 
crete ones (black dots): the plots on the top correspond to the variables q\ — Q (on the 
right) and pi (on the left): the values of the initial conditions and of the parameters Ao 
and \x used are the same as for figure ([1]). For the two plots on the bottom we changed 
only the value of \i from 5 to 2: as explained in Section (J2]) and as can be seen from 
the figures, to a smaller value of the parameter /i it corresponds a smaller value of the 
time step size. 
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